Vibrating-coil magnetometry of the spin liquid properties of Tb 2 Ti 2 07 



S. Legl, 1 C. Krey, 1 S. R. Dunsiger, 1 H.A. Dabkowska, 2 J. A. Rodriguez, 3 ' 4 G.M. Luke, 3 and C. Pfleiderer 1 

1 Technische Universitat Miinchen, Physik- Department E21, D-85748 Garching, Germany 
2 Brockhouse Institute for Materials Research, McMaster University, Hamilton, Ontario L8S l^Ml, Canada 
Department of Physics and Astronomy, McMaster University, Hamilton, Ontario L8S 4-M1, Canada 
* Laboratory for Muon Spin Spectroscopy, Paul Scherrer Institut, 5232 Villigen, Switzerland 

(Dated: March 8, 2013) 

We have explored the spin liquid state in Tb2Ti207 with vibrating coil magnetometry down to 
~ 0.04 K under magnetic fields up to 5 T. We observe magnetic history dependence below T* ~ 0.2 K 
reminiscent of the classical spin ice systems Ho2Ti207 and Dy2Ti2C>7. The magnetic phase diagram 
inferred from the magnetization is essentially isotropic, without evidence of magnetization plateaux 
as anticipated for so-called quantum spin ice, predicted theoretically for [111] when quantum fluctua- 
tions renormalize the interactions. Instead, the magnetization for T<T* agrees semi-quantitatively 
with the predictions of "all- in/all-out" (AIAO) antiferromagnetism. Taken together this suggests 
that the spin liquid state in Tb2Ti207 is akin to an incipient AIAO-antiferromagnet. 

PACS numbers: 75.30.Kz; 75.60.Ej; 75.40.Cx; 75.40.Gb 



Pyrochlore oxides, A2B2O7, in which rare earth mag- 
netic moments are located on the A-site of a three- 
dimensional network of corner sharing tetrahedra are 
model systems of geometric frustration pQ. The con- 
sequences of such geometric frustration are intimately 
connected with the strength and the nature of the mag- 
netic anisotropy at the rare earth site. For instance, 
in Ho2Ti2C>7 and Dy2Ti207 a strong easy-axis (Ising)- 
anisotropy along the local [111] axis in the unit cell, to- 
gether with net ferromagnetic interactions, are the most 
important preconditions for the emergence of the highly 
celebrated spin ice behavior [2HS]. 

A major unresolved question in geometric frustration 
concerns the fate of the spin ice state, when the strength 
of the local Ising anisotropy is reduced. This may boost 
the relative importance of quantum fluctuations. In fact, 
an exciting theoretical proposal states that quantum fluc- 
tuations may renormalise the exchange interactions of 
an unfrustrated (111) antiferromagnetic Ising antiferro- 
magnet, making them effectively ferromagnetic. The as- 
sociated novel state is referred to as quantum spin ice 
(QSI) [6-8], and may be viewed as a spin liquid in which 
thermal longitudinal spin fluctuations which break the 
spin ice rules, as well as thermal and quantum fluctua- 
tions transverse to the local (111) directions are relevant. 
Compelling evidence of a QSI would be magnetization 
plateaux like those observed in Hx^T^Cv and Dy2Ti2C>7, 
for a magnetic field strictly along a global (111) axis 0- 
IT2"] . In turn an anisotropy of the magnetic phase diagram 
would be expected reminiscent of classical spin ice. 

An ideal model system to study whether any of the 
classical spin ice properties survive under reduced local 
magnetic anisotropy is Tb2Ti2C>7. At high temperatures 
Tb2Ti2C>7 exhibits a Curie-Weiss susceptibility with a 
large effective moment fj, c g = 9.6/ieTb -1 and a nega- 
tive Curie- Weifi temperature ©cw characteristic of an- 
tiferromagnetic interactions [131 03] ■ As the crystal elec- 



tric field (CEF) of the Tb 3+ ion leads to a ground state 
doublet and an energy gap of ~ 18 K to the first ex- 
cited state [Mj [15] , antiferromagnetic order is expected 
around ~ IK [To] . Surprisingly however, /^SR [T51 [T7] . 
the ac susceptibility [T^J [TS], and Neutron Spin Echo 
(NSE) [JSJ [Jj5] established strong spin dynamics down 
to 20 mK without long-range magnetic order. 

The origin and nature of the lack of long-range mag- 
netic order in Tb2Ti207 represents a major puzzle in 
geometrically frustrated magnetism. It must however be 
a sensitive function of the low lying energy levels. X-ray 
diffraction [20] and a related transition at ~0.15K [2T] 
suggest a cooperative Jahn- Teller distortion. A possible 
splitting of the doublet to yield a singlet ground state 
has been inferred from specific heat data [32], though 
this is at odds with high energy resolution neutron- 
scattering [23.. The large electronic-nuclear hyperfine 
coupling of Tb suggests that nuclear degrees of freedom 
may also be important. It is unresolved if evidence for 
magnetic glassiness in the milli-kelvin regime is intrin- 
sic [24] or due to defects [HI [25]. Moreover, for mag- 
netic fields applied along (110) the specific heat [26] and 
neutron scattering [27] suggest field-induced order with 
spin-ice like properties |28H30j . while the spin order for 
field along (111) is undetermined [3TJ. Magneto-elastic 
effects are also known to be very pronounced [20J [32j [33] , 
consistent with the observation of pressure induced anti- 
ferromagnetic order |14j . 

Calculations exploring the role of the interaction 
strength and CEF splitting suggest that Tb2Ti 2 07 is at 
the border between a QSI and "All In All Out" (AIAO) 
antiferromagnetism [5HS]- However, because the zero- 
frequency QSI correlations may be masked by strong 
fluctuations and difficult to detect experimentally, it has 
been emphasized that magnetization measurements at 
mK temperatures are ideal to identify a QSI, since they 
probe the zero frequency (and zero wave vector) response. 
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The possible existence of QSI in Tb2Ti207 has been 
addressed by Baker et al. [34] , who infer the existence of 
magnetization plateaux from ac susceptibility data down 
to 25 mK, as well as /iSR measurements which exhibit 
peaks and kinks in the spin lattice relaxation rate consis- 
tent with the boundaries of the magnetization plateaux. 
In another study, l'Hotel et al. (35] recently report ex- 
traction magnetometry and ac susceptibility down to 
80 mK claiming the absence of magnetization plateaux. 
However, data reported in the former study probe the 
response at finite frequency, while the latter study does 
not extend to low enough temperatures to be conclusive. 
Moreover, extraction magnetometry suffers from the risk 
of tiny sample vibrations in the magnetic field, thereby 
changing the field history of the sample. Finally, both 
studies addressed a (111) direction only, not reporting 
the properties of other directions as a control experiment. 
Taken together, the nature of the spin liquid state in 
Tb2Ti2C>7 and the proposal of QSI are hence unresolved. 

In this Letter we address, to the best of our knowl- 
edge for the first time, the properties of Tb2Ti 2 07 and 
the possible existence of QSI down to sufficiently low 
temperatures under conditions avoiding parasitic mod- 
ifications of the field history. The magnetization of a 
Tb 2 Ti 2 07 single-crystal was measured at TUM using a 
bespoke vibrating-coil magnetometer (VCM) for temper- 
atures down to ~ 0.04 K and magnetic fields up to 5 T 
|36|, I37j . Our VCM represents a new development previ- 
ously thought to be prohibitively difficult. As the main 
advantage of the VCM, the sample does not move during 
measurements (the detection system is completely decou- 
pled from the dilution refrigerator). The coil-set was op- 
erated at 36.5 Hz and the sample temperature was mon- 
itored with several RuC>2 and Speer sensors. Great care 
was taken to ensure good thermal anchoring (for details 
see [HIST]). 

Magnetic field sweeps at up to IT were recorded in 
a step mode, where the field was kept constant while 
recording the magnetization. Magnetization data up to 
5 T were recorded while sweeping the magnetic field con- 
tinuously at 15 mT min -1 . For studies of a (100) and 
(110) direction the sample was mounted on a small cop- 
per wedge and the orientation confirmed with Laue x-ray 
diffraction. The empty sample holder was measured care- 
fully and the background signal subtracted. The VCM 
data recorded below 1.5 K was calibrated by means of 
a Ni standard measured in the same parameter range. 
A comparison to the extrapolated magnetization of the 
Tb2Ti207 sample measured at 1.8 K in a Quantum De- 
sign MPMS system enabled the background contribution 
to be isolated. High temperature data were recorded in 
an Oxford Instruments vibrating sample magnetometer. 

For our study a single crystal was grown at McMas- 
ter University by optical float-zoning. The feed and seed 
rods were prepared by annealing high-purity Tb407 and 
Ti02 in air. The sample was then float-zoned in a high- 
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FIG. 1: (Colour online) Temperature dependence of the mag- 
netization of Tb2Ti207 in a small applied magnetic field of 
10 mT. Below T* ~ 200 mK a distinct difference between data 
recorded under zero-field cooling (zfc) and field-cooling (fc) 
emerges, characteristic of a glassy magnetic state (all data 
were recorded while field-heating (fh) to minimise systematic 
errors). Curves are shifted for clarity as described in the text. 

purity Ar atmosphere of 4 bar with a rate of 7mmhr -1 
|38) . Powder x-ray diffraction of a small piece of the 
single crystal confirmed that the ingot was phase-pure 
with the correct crystal structure. A single-crystalline 
disc (7 x 4.1 x 0.8 mm 3 ) was cut from the ingot for the 
magnetization measurements. The disc was oriented per- 
pendicular to a (111) direction within ~ 1°. The mag- 
netic properties of this crystal as recorded above ~ 1.5 K 
were in excellent agreement with the literature. More- 
over, none of our data displayed anomalies in dM/dT and 
thus the magneto-caloric effect, probing spurious features 
in the specific heat seen in low-quality Tb2Ti2 07 |2"2"] . 

Demagnetizing fields were corrected for the [111] ori- 
entation, approximating our sample as an ellipsoid with 
a demagnetising factor N — 0.34. As we will show, the 
conclusions drawn from our data for [100] and [110] do 
not depend on the correction of demagnetizing fields, so 
only raw data as a function of applied field are shown. 
The internal fields were estimated to deviate by a few 
degrees from the [100] and [110] directions, because the 
plane of the disc was tilted with respect to the applied 
field. In addition, the demagnetization fields in the tilted 
sample may have been slightly inhomogeneous. 

Figure [T] illustrates the temperature dependence of the 
magnetization in an applied field of 10 mT, where data 
for [100] and [110] have been shifted by 0.03^ B Tb" 1 and 
0.06/iBTb -1 , respectively for clarity (the uncertainty in 
the constant background is ~ ±0.015//b Tb _1 ). Prior 
to recording these data the superconducting magnet was 
demagnetized at ~ 2.3 K to remove any parasitic rema- 
nent fields. After zero-field cooling the field was increased 
at a rate of lmTmin -1 to the setpoint of 10 mT and 
data recorded while heating the sample continuously with 
5mKmin _1 up to ~ 1.5 K (zfc-fh). Following this the 
sample was cooled down with the field unchanged and 
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data recorded while heating at the same rate (fc-fh). 

With decreasing temperature the magnetization in- 
creases with a positive curvature consistent with the 
paramagnetic properties at high temperatures. In all 
field directions, the curves display a cusp in the zfc- 
fh and fc-fh data. We find slightly different values of 
T* ~ 0.246 K, - 0.212 K and - 0.204 K for [110], [100] 
and [111], respectively not reported previously. For 
T > T* the zfc-fh and fc-fh data agree essentially (the 
nature of the tiny difference for [110] is most likely asso- 
ciated with a small drift of the detection system). 

The shape of the cusp, the absolute difference of zfc-fh 
and fc-fh data and the qualitative temperature depen- 
dence of the data provide strong evidence of the emer- 
gence of intrinsic magnetic glassiness below T* which is 
essentially isotropic. Note that magnetic ordering coex- 
isting with persistent low temperature fluctuations has 
been reported in a variety of geometrically frustrated sys- 
tems [39l 00] . It is therefore perfectly consistent with the 
persistent muon spin relaxation observed in in Tb2Ti207 
using /xSR [T31 |T7] . The absolute size of the magneti- 
zation for T > T* , being largest for [110] and smallest 
for [111], suggests that the tiny variation in T* and the 
difference between fc-fh and zfc-fh, which are also largest 
for [110] and smallest for [111], originate from the mag- 
netic anisotropy. However, in comparison to simple cubic 
systems the anisotropy is unusual, as [110] cannot be a 
soft direction. 

Presented in Figs.[2](A), (B), (C) and (D) are magneti- 
zation data in the range from ~ 0.043 K to ~ 50 K under 
magnetic fields up to 5 T for [110], [100] and [111], re- 
spectively. Data below 1 K were recorded after zero-field 
cooling. For the field and temperature scale shown in 
Fig.[2]field-cooled data are identical. Note that panel (D) 
shows the data of panel (C) as a function of calculated 
internal field. Apart from a shift of the characteristic fea- 
tures in the magnetization towards lower fields, demag- 
netizing fields do not affect the conclusions of our study. 
For the temperature and field scale shown here data are 
qualitatively rather similar. The linear magnetic field 
dependence at ~ 50 K becomes highly non-linear at the 
lowest temperatures studied. It remains unsaturated at 
5 T, even though the magnetization reaches a large value 
between 5 and 6 /ib Tb _1 . 

On the scale shown in Fig. [2] we find no evidence for a 
magnetization plateau, predicted for QSI. This contrasts 
with classical spin ice systems, where the magnetization 
plateaux are a prominent feature strictly along (111) on 
the same scale. Moreover, we observe some fine-structure 
in dM/d/Xo-ff not reported before (Figs.[2](E), (F), (G) 
and (H)). Namely, for the lowest temperature the cal- 
culated first derivatives allow us to define two cross-over 
scales /in Hi and fio H 2 . The field (j,q Hi marks the end of 
the initial rise of the magnetization associated with the 
initial drop of dM/d/ioH , while /in Hi > /in Hi marks 
a faint drop of dM/dfioH. Thus the magnetization for 
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FIG. 2: (Color online) Magnetization of Tb 2 Ti 2 7 for tem- 
peratures in the range ~ 0.043 K to ~ 50 K under magnetic 
fields up to 5 T. Data below 1 K were recorded after zfc. Pan- 
els (A), (B) and (C) display data as a function of the applied 
field along [110], [100] and [111]. Panel (D) shows the data 
of panel (C) as a function of internal field. Panels (E), (F), 
(G) and (H) show the numerical derivative of the experimental 
data recorded at the lowest temperatures in order to illustrate 
the definition of the characteristic fields /in Hi and /in H2 ■ 



/in Hi < B < /in H2 cannot correspond to a plateau (or 
remnants thereof), since the slope decreases when ex- 
ceeding fiQ H 2 without a point of inflection at /in H 2 . In 
addition, /i Hi and /in H 2 exist for all directions and are 
not specific for [111]. 

Based on the isotropic behaviour observed, we specu- 
late that microscopic features of /in Hi and /in H 2 may 
have been seen in time-of-flight neutron scattering for 
field parallel (110), where the magnetic diffuse scattering, 
condenses into a new Bragg peak around /in Hi consis- 
tent with a polarized paramagnet |27j . A magnetically 
ordered phase, which supports spin wave excitations is 
induced around /in H 2 , consistent with cross-over phe- 
nomena around in the specific heat and ac susceptibil- 
ity [27] . Above /i H 2 magnetization plateaux character- 
istic of spin ice are therefore clearly no longer expected. 

Shown in Fig.[3](A) are typical zero-field cooled mag- 
netization data for the [111] direction, where we find no 
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FIG. 3: (Color online) Experimental and theoretical low-field 
magnetization of Tb2Ti207 for the [111] axis in the parameter 
range of the predicted magnetization plateau of quantum spin 
ice. Only zero-field cooled data are shown. (A) Magnetic field 
dependence of the magnetization of Tb2Ti207 in small fields 
at various temperatures. (B) Numerical derivative of the data 
shown in panel (A). (C) Theoretically predicted magnetiza- 
tion for a quantum spin ice state in Tb2Ti207 as reported in 
Ref. [8] . (D) Derivative of the theoretical data shown in panel 
(C). 



evidence of magnetization plateaux either. Qualitatively 
field-cooled data are the same as shown in the supplemen- 
tary material. The absence of magnetization plateaux 
is most evident in dM/dB- m t calculated from the data 
shown in Fig.[3](B). This plot does not display a point 
of inflection of M(B mt ). Instead, dM/dB mt has a broad 
maximum only. Data for [110] and [100] shown in the 
supplement are similar to [111]. Our data are thereby 
consistent with Ref. [35] , but extend a factor of two lower 
in temperature, well into the proposed QSI regime. 

For comparison with the experimental data we repro- 
duce in Fig.[3](C) and (D) theoretical calculations of the 
magnetization and their first derivatives at 20 mK for 
J = 0.167 ("quantum spin ice" (QSI)) and J = 0.2 
("all-in/ all-out" (AIAO)) 0. The most important qual- 
itative difference between the QSI and AIAO concerns 
the marked change in the slope of dM/d^H at low 
field (< 0.05 T) from negative to positive for QSI and 
AIAO structures, respectively. Notice that experimen- 
tally AM/ dB[ nt at 43 and 70 mK has a positive slope at 
low field and a maximum at 0.03 T. This observation is 
independent of the field history as described above. 

However, the calculated AIAO spin state is based on 
a crystal field splitting 1/ A which is reduced by a factor 
of two as compared with experiment and an exchange 
coupling, J = 0.2, that is larger than that inferred ex- 
perimentally from the AIAO spin structure in large mag- 
netic fields [14]. Unfortunately we cannot offer an ex- 
planation for the discrepancy of the crystal field split- 
ting. As for the exchange coupling it is important to 
note that the experimental value was determined in large 
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FIG. 4: (Color online) Magnetic phase diagram of Tt>2Ti207 
for [110], [100] and [111] as inferred from the magnetization. 
The cross-over fields /in Hi and /in H2 are larger the magneti- 
cally softer the direction. This may serve as a testing ground 
for the validity of theoretical models. Open symbols in panel 
(C) correspond to the phase diagram after correction of de- 
magnetising effects. 



magnetic fields and may therefore differ from low fields 
|14j . Moreover, as shown in Ref. 7] for J cx = 1/6 K the 
predicted ground state depends on the model and its im- 
plementation [7] , where the dipolar spin ice model and a 
cubic unit cell model with Ewald summed dipole-dipole 
interactions and crystal electric fields predict AIAO an- 
tifcrromagnetism and long-range spin ice with ordering 
wave vector q = (000), LRSLjoO) respectively. Note how- 
ever that the phase diagram is sensitive to the model 
used [H]. I n particular, the model of Refs. [6-8] com- 
pletely ignores the possibility of anisotropic exchange and 
higher multipolar couplings. Because Tb2Ti 2 07 does not 
develop spontaneous long-range magnetic order, it seems 
natural to conclude that the spin liquid in Tb2Ti20y is 
dominated by strong fluctuations at the border of AIAO- 
antiferromagnetism. This may be referred to as incipient 
AIAO-antiferromagnetism. 

The nearly isotropic magnetic phase diagrams inferred 
from our magnetization data for [110], [110] and [111] 
shown in Fig. [4] corroborate incipient AIAO antiferro- 
magnetism. Note that the open symbols in panel (C) 
show the phase diagram when correcting demagnetising 
effects. The expected corrections in panels (A) and (B) 
are smaller than for panel (C). As a function of magnetic 
field we find two cross-over scales Hi and /io H 2 , where 
neutron scattering suggests different changes of the un- 
derlying microscopic properties [2"T] . 

In conclusion, we find no evidence of magnetization 
plateaux in Tb2Ti20y expected of QSI and fluctuation- 
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induced ferromagnetic interactions. Instead, our data 
are in semi-quantitative agreement with the theoretical 
predictions of AIAO-antiferromagnetism, suggesting that 
the spin liquid state in Tb 2 Ti 2 07 may be viewed as an in- 
cipient AIAO antifcrromagnet. The small remaining ori- 
entational dependence in T* , fi Hi and H 2 reported 
in this Letter, which corresponds with the magnitude of 
the magnetization and thus the magnetic anisotropy, no 
doubt provides an important clue in future studies as 
to what inhibits the formation of long-range AIAO anti- 
ferromagnetism in Tb2Ti20y and drives the spin liquid 
state microscopically. 
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SUPPLEMENTARY MATERIAL FOR: 
VIBRATING-COIL MAGNETOMETRY OF THE 
SPIN LIQUID PROPERTIES OF TB 2 TI 2 7 

Shown in Fig. [5] are low temperature magnetization 
data of our Tb2Ti2C>7 sample. For clarity only two tem- 
peratures are shown, one well below and the other well 
above T*. For these measurements the VCM and sam- 
ple were at first demagnetised at 2.3 K and subsequently 
cooled to the desired temperature. Data were recorded in 
a sequence of three field sweeps at constant temperature: 
(i) from fi H = to +1 T, (ii) from +1 T to -1 T, and 
(iii) from —1 T to +1 T. Thus, the first sweep represents 
the zero-field cooled magnetization and the second and 
third sweep the field-ccoled magnetization. 

Data for field along [110] and [100] are shown in 
Fig.[5](A) and (B) as a function of applied magnetic field 
without correction for demagnetising effects. Data for 
[111] are shown as a function of internal field in Fig. [51(C) 



taking into account demagnetising fields. The fine struc- 
ture of the slope of the magnetization is best revealed in 
the calculated first derivative of the magnetization. At 
the lowest temperatures measured the derivative displays 
a small maximum, regardless of sample orientation and 
field history in contrast to the prediction of a QSI. 




FIG. 5: (Color online). Typical magnetization data of our 
Tfc>2Ti2 07 data comparing zero-field cooled with field cooled 
behaviour. Qualitatively no difference is observed regardless 
of orientation. 



